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ABSTRACT i ^  Mode-I  crack  growth  in  an  elastic  perf ectly-plastic  material 
under  conditions  of  generalized  plane  stress  has  been  investigated.  In 
the  plastic  loading  zone,  near  the  plane  of  the  crack,  the  stresses  and 
strains  have  been  expanded  in  powers  of  the  distance,  y,  to  the  crack 
line.  Substitution  of  the  expansions  in  the  equilibrium  equations,  the 
yield  condition  and  the  constitutive  equations  yields  a  system  of  simple 
ordinary  differential  equations  for  the  coefficients  of  the  expansions. 
This  system  is  solvable  if  it  is  assumed  that  the  cleavage  stress  is 
uniform  on  the  crack  line.  By  matching  the  relevant  stress  components 
and  particle  velocities  to  the  dominant  terms  of  appropriate  elastic 
fields  at  the  elastic-plastic  boundary,  a  complete  solution  has  been 
obtained  for  'ey  in  the  plane  of  the  crack.  The  solution  depends  on 
crack-line  position  and  time,  and  applies  from  the  propagating  crack 
tip  up  to  the  moving  elastic-plastic  boundary.  Numerical  results  are 
presented  for  the  edge  crack  geometry.  ^ — 

£^i!^  f,J'9  'i 
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Introduction 

Quasi-static  fields  of  stress  and  deformation  near  the  tip  of  a 
growing  crack  in  an  elastic  perf ectly-plastic  material,  have  been 
discussed  in  considerable  detail  by  Rice  [1] .  Analytical  expressions 
for  the  near-tip  fields  reveal  the  asymptotic  structure  of  the  fields  in 
the  immediate  vicinity  of  a  moving  crack  tip.  They  contain,  however, 
functions  that  can  generally  be  obtained  only  by  supplemental  numerical 
procedures.  The  one  complete  analytical  solution,  which  is  for  the 
case  of  anti-plane  strain,  has  been  given  by  Rice  [2], 

In  a  recent  paper  Achenbach  and  Dunayevsky  [3]  considered  the  case 
of  Mode-I  crack  growth  under  plane  stress  conditions  in  an  elastic 
perfectly-plastic  material.  In  Ref. [3]  it  was  assumed  that  the  stress 
components  for  a  centered  fan  field,  which  were  discussed  by  Hutchinson  [4], 
and  which  satisfy  the  yield  condition  and  the  equilibrium  equations,  are 
valid  up  to  the  elastic-plastic  boundary  (at  least  near  the  crack  line) . 

The  analytical  approach  of  Ref. [3]  then  employs  expansions  of  the  particle 
velocities  in  powers  of  y  (the  distance  from  the  plane  of  the  crack) ,  to 
obtain  ordinary  differential  equations  with  respect  to  x  for  the  co¬ 
efficients  in  the  expansions.  Functions  of  time  that  enter  in  integrating 
these  equations  were  determined  by  matching  the  fields  in  the  plastic 
loading  zone  to  the  dominant  terms  of  suitable  elastic  fields  at  the 


elastic-plastic  boundary. 


In  the  present  paper  we  reconsider  the  results  of  Ref. [3]. 

Instead  of  making  an  a-priori  assumption  on  the  stress  field,  we  also 
write  expansions  for  the  stresses  near  the  crack  line.  Substitution 
of  these  expansions  into  the  equilibrium  equations  and  the  yield 
condition  produces  a  simple  system  of  equations  for  the  coefficients. 
Unfortunately  the  system  is  not  closed,  and  additional  information  is 
required.  The  structure  of  the  equations  suggests  that  the  cleavage 
stress  on  the  crack  line  is  uniform.  With  that  assumption,  the  system 
of  equations  can  be  solved  in  a  simple  manner.  The  resulting  expressions 
for  the  stresses  are  consistent  with  the  centered  fan  field  employed 
in  Ref .[3] . 

To  illustrate  the  analytical  results  the  plastic  strains  just 
ahead  of  a  moving  crack  tip  have  been  computed  for  an  edge  crack 
geometry.  The  results  are  particularly  suited  for  use  in  conjunction 
with  a  critical  strain  criterion. 


Governing  Equations 

The  geometry  that  is  being  considered  in  this  paper  is  shown  in 
Fig.  1.  The  x^-axis  of  a  stationary  coordinate  system  is  parallel  to 
the  crack  front,  and  x^  points  in  the  direction  of  crack  growth.  The 
position  of  the  crack  tip  is  defined  by  x^  =  a(t).  A  moving  coordinate 
system,  x,y,z  is  centered  at  the  crack  tip,  with  its  axes  parallel  to 
the  X^»X2  anc*  x3  axes*  Relative  to  the  moving  coordinate  system  we  also 
define  polar  coordinates  r,0,  with  0*0  coinciding  with  the  positive 
x  direction. 

In  the  moving  coordinate  system  the  equilibrium  equations  are 
3o  3t  3t  3o 


3x 


jnc 


3y 


-.*z  +  —L 

3x  3y 


(la.b) 


We  consider  a  state  of  generalized  plane  stress,  hence  o  ,  o  and  o 

z  xz  yz 

vanish  identically .  The  Huber-Mises  yield  criterion  may  then  be 
written 


o2  +  a2  -  a  o  +  3t2 


3kz 


x  y  x  y  xy 
where  k  is  the  yield  stress  in  pure  shear.  The  strain  rates  are 


(2) 


•  3u  .  3v  .  l,3u  3v. 

Ex  3x  *  Ey  3y  ’  Exy  2  3y  3x 

In  the  moving  coordinate  system  the  material  time  derivative  is 


(3a,b,c) 


where  a  =  da/dt  is  the  speed  of  the  crack  tip.  The  strain  rates 
are  related  to  the  stresses  and  stress  rates  by 


(4) 


(5) 


=  T’Cd  -  vj  )  +  i  A  (2a  -  a  ) 

9x  E  x  y  3  x  y 


=  i(<j  -  vo  )  +  -i-  a (2a  -  a  ) 

9y  E  y  x  3  y  x 


(6) 


1 ,9u  9v. 

2{dy  9x ' 


1+v 


t  +  At 
xy  xy 


(7) 


where  E  and  v  are  Young's  modulus  and  Poisson's  ratio,  respectively 
and  A  is  a  positive  function  of  time  and  the  spatial  coordinates. 


Solution  along  the  Crack  Line 

In  this  paper  we  are  irterested  in  solutions  along  the  crack  line 
y  =  0,  0  <  x  x^  where  x  =  x^  defines  the  elastic-plastic  boundary. 
Such  solutions  can  be  obtained  by  considering  expansions  with  respect 
to  y  in  the  region  y/x  <<  1: 


°x  =  +  P2(x,t)yz  +  p4(x,tjy4  +  ... 

oy  “  qQ(x,t)  +  q2(x,t)y2  +  q4(x,t)y**+  ... 

=■  s1(x,t)y  +  s3(x,t)y3  +  ... 

u  =  uQ(x,t)  +  u2(x,t)yz  +  ... 
v  -  v](x,t)y  +  (x, t)y3  +  ... 

A  *  AQ(x,t)  +  A2(x,t)yz  +  ... 

Here  we  have  taken  into  account  that  u  and  A  are  symmetric 

,  v  are  antisymmetric.  Substition  of 


(8) 

(9) 

(10) 

(11) 

(12) 

(13) 


with  respect  to  y  =  0,  while  a 
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(8)-(10)  into  (la,b)  and  collecting  terms  of  the  same  order  in  y  yields 


+  s.  =0 
3x  1 


+  3s.  =  o 

3x  3 


ST  +  2q2  =  0 


ir  +  4q4 =  0 


Substitution  of  (8) -(10)  into  the  yield  condition  (2)  yields  by  the  same 
procedure 


p2  +  q2  -  p  q  -  3k2 
*o  ^o  o  o 


(2P0-q0)P2  +  (2Vpo)q2  +  3sl  =  0 


P2  +  (2po-qo)p4  +  q2  +  (2qo’Po)q4  ’  P2q2  +  6slS3  =  ° 
In  the  same  manner  we  obtain  by  using  (8) -(13)  in  (5) -(7) 


3u  .  3p  3q  .  3p  3q  A 

_ o  1,  o  oN  a,  o  ,,  Os  .  o 

3x  =  El3t  3t  ^  E  3x  3x  ^  +  3  (2po_CV 


^2  ,  1.^2 

3x  =  E'3t 


U  a  3p2  3q2  Ao 

/ _ A  . . _ A\  i  _ o  t  • i 


(3 r  - v  ir>“  Pit  - v  *r>  +  (Zp2  -  V  +  -7  (2W  <22> 


.  3q  3p  .  3q  3p  A 

•  1  /  o  ..  O'  a ,  o  o .  .  o  N 

V1  “  E(3t  "  3t  )  "  E  3x  "  V  3x  ;  3  (2qo“Po) 


Br 


pi*J 


3*3  -  flf  -  V  if)-  f(sf  '  V  sf )  +  IAo(2W  +  ?V2W  (24) 


i  3v,  .  3s,  3s, 

.  ,1  1  1+v  , 1  .  1.  ,  . 

uo  +  T  "o =  r  vTI - a  T )  +  A  S. 

2  2  3x  E  3t  3x  o  1 


where  (4)  has  also  been  used, 


At  this  stage  we  have  14  unknowns  and  12  equations.  It  turns  out 
that  this  system  of  equations  can  be  solved  provided  that  one  assumption 
is  made.  We  assume  that  a  is  constant  on  the  crack  line  in  the  plastic 

y 

loading  zone.  Hence,  by  the  use  of  Eq.(9): 

q  =  constant  (26) 

o 

It  will  be  shown  later  that  the  assumption  is  consistent  with  the 
stress  components  for  a  centered  fan  field  which  have  been  considered  by 
Hutchinson  [4],  and  which  are: 


a  =  kcos30,  a  =  k(2cos38  +  3sin29cos0),  a  =  -ksin30 
x  y  xy 


(27a, b.c) 


If  q  =  constant,  it  follows  from  Eq.(18)  that  p  =  constant, 
o  o 

On  the  basis  of  p^  =  constant,  it  follows  from  Eq.(14)  that  s^  =  0, 
and  subsequently  from  Eq.(16)  that  q^  =  0.  Substitution  of  these 
results  in  Eq.(19)  yields  q^  =  2pQ,  anc*  Eq.(18)  then  gives  pq  =  k  and 
qQ  =  2k.  Equation  (20)  subsequently  yields  q^  ■  -p2/3k.  Substitution 
of  the  latter  result  in  Eq.(17)  and  then  in  Eq.(15)  yields 


The  solution  to  Eq.(28)  which  satisfies  the  condition  that  the  stresses 


are  multivalued  at  the  crack  tip  is 


In  summary,  it  has  been  shown  that  the  equilibrium  equations  and  the 


yield  condition  are  satisfied  by 


_  1(L\  i  +  nA 


ax  =  k[l  -  f  (J>  ]  +  0(J) 


a  -  2k  +  0(z) 
y  x 


T  =-k(I)3+0(i)5 

xy  x  x 


It  is  noted  that  Eqs . (30)- (31)  are  indeed  expansions  with  respect  to 
0  «  y/x  of  Eqs. (27a, b,c) . 

In  the  next  step  we  substitute  Eqs . (30)-(32)  into  (21)-(25)  and 
collect  terms  of  the  same  order.  The  result  is 


fV  „  ,  k  i  i  • 

3x  ’  3x  E  jc3"  xT  o 


v  -  A  k 
1  o 


3v. 


2"2  +1T 


By  combining  03)  .04)  and  05)  we  obtain 


i  8  v.  v  .  . 

1  _ 1  1  _  -  k  a 

2  3x  J  E? 


The  general  solution  to  Eq.(36)  is 


v .  -  |  {-2  ^  4n(5-)  +  +  C (t)x2 }  , 


(37) 


8 


where  x^  defines  the  x-coordinate  of  the  elastic -plastic  boundary 
on  the  crack  line.  The  functions  B(t)  and  C(t;  can  be  obtained 
from  continuity  conditions  at  the  elastic-plastic  boundary. 

For  small  values  of  0  (i.e.,  y/x  <<  1)  the  field  in  the  plastic 
loading  zone  will  be  matched  at  the  elastic-plastic  boundary  to  the 
dominant  terms  of  a  corresponding  elastic  field.  For  the  elastic  field 
we  do,  however,  not  take  the  field  for  a  crack,  but  rather  that  for  a 
notch  with  as  radius  of  curvature  at  its  tip.  In  polar  coordinates 
R,y,  the  appropriate  Mode-I  stress  fields  are  given  by  Creager  and  Paris  [  5] 


°x  =  (^r)  hicoa¥  S±nf*]  "  2R  C°S¥i  ' 

°y  ‘  (^r)  KIlC°4>  I1+Sl4>  »in|*1  +  2R  C°Sf»'1  1 

Txy  *  (lk)  VSi41'  c°4,,':OS2*  •  2R  1 

Note  that  the  tip  of  the  notch,  which  is  not  the  tip  of  the  crack 
nor  the  elastic-plastic  boundary,  is  a  distance  from  the  origin  E , 
as  shown  In  Fig.  1.  The  center  of  the  elastic  field  E,  whose  position 

is  defined  by  x^  =  e(t),  y^  =  0,  is  located  in  between  the  crack  tip 

and  the  elastic-plastic  boundary  defined  by  x  =  x^tt).  For  generalized 
plane  stress,  the  displacements  corresponding  to  (3b)-(40.)  are 

h 

u  =  iu  ki^cos2  ~  l+2sin2-|^]  +  cos^}  < 


9 


notch  rather  than  for  a  crack.  For  an  elastic  crack-tip  field  the 
conditions  of  reaching  the  yield  condition  at  the  elastic-plastic 
boundary  would  conflict  with  the  condition  of  continuity  of  o^,  as 
can  be  checked  by  setting  p  =  0  in  (43)  and  (44). 


10 


Further  details  of  the  matching  procedure  have  been  given  by 
Achenbach  and Dunayevsky  [3  ].  The  relevant  results  are 


x„  =  (l-y)x  ,  and  thus  R  =  yx  , 
E  p  P  P 


(47, a, b) 


where  y  =  1/ Jl . 


B(t)  =  B.a(t)  +  B0x  (t)  (48) 

1  2  p 


C(t)  =  [(^(t)  +  C2xp(t)]/[xp(t)]3, 

(49) 

where 

Bi  =  33  V  [K+5  +4y(k+1)J(e/v)  -  | 

(50) 

B2  =  33  y7  +  2y(k+1)]  (E/v) 

(51) 

C.  "  77T-  1  "(<+5)  +  2y  (<+l)  ]  (E/p)  +| 

1  32  y  3 

(52) 

C2  =  32  y*  f_(K+5)  +  4y (k+1)](E/v) 

(53) 

The  Strain  on  the  Crack  Line 

In  the  plane  of  the  crack  we  have  =  v^.  At 

the  elastic  plastic 

boundary  (42)  yields  for  small  y 

3k.  1.  2-v  . 

^y^PB  8  p  3^  '  E 

where  (46),  (47)  and  8  ~  y/x  have  been  used.  In  the  stationary 

coordinate  system,  (37)  can  now  be  integrated  to  yield  the  total 

strain  for  t  >  t  as 
~  P 


(54) 


where  v^(x^,s)  is  obtained  from  (37)  by  using  the  relation 

x  =  x,  -  a(t).  In  (55),  t  is  the  time  that  the  elastic-plastic 
1  P 

boundary  arrives  at  position  x^.  Thus,  for  the  propagating  crack  tip, 
t  follows  from  the  equation 

a(t  )  +  x  (t  )  =  x  ,  (56) 

P  P  P  1 

Here  x  (t  )  is  obtained  from  the  stress  intensity  factor  by  using 
P  P 

(46)  and  (47b): 

xP(t)  =  T»  tVt)/k]2  (57) 

Equation  (37)  is  also  valid  for  a  stationary  crack.  By  setting 

a  =  0  we  obtain  from  (37)  for  t  >  t 

—  p 


vfC(A,t)  =  ^  [B-/A  +  C_A2/x3(t)]  x  (t)  ,  (58) 

1  E  i.  l  p  p 

where  (48)  and  (49)  have  been  used  and  A  is  a  distance  ahead 
of  the  crack  tip.  The  arrival  time  of  the  elastic-plastic  boundary 
follows  from  (56)  as 


a  +  x  (t  )  =  x- 
o  p  p  1 


Let  us  consider  the  case  that  loading  starts  at  time  t  =  0, 


but  that  the  crack  tip  does  not  start  to  propagate  until  time  t  *  t  . 

s 

For  a  position  x.  which  is  inside  the  plastic  zone  at  time  t  =  t  ,  we 
-L  s 


12 


Ey(xl,l:^“  (ey)PB  +  /  v®C(xltt)ds  +  /  v^.s^s 


Here  t  follows  from  (59),  and  t  Is  the  time  that  the  crack  tip  arrives 
P  e 

at  a  small  distance  A  from  the  position  x^  that  is  being  observed.  We 


x  -  A  =  a(t  ) 
l  e 

For  a  position  which  is  outside  the  plastic  zone  at  time  t  =  tg. 


X1  >  ao  +  xp(ts)  *  (62 

the  expression  (55)  holds,  where  t  is  now  defined  by  (56). 

P 

Equations  (55)  and  (60)  can  be  manipulated  to  yield  the  singular 
parts  of  the  strain  at  x^  =  a(t),  plus  a  bounded  integral.  The 
result  can  be  found  in  Ref. [3].  In  the  present  paper  the  integrals 
(55)  and  (60)  have  been  evaluated  numerically. 

The  result  simplifies  considerably  for  the  case  that  all  fields 
are  assumed  to  be  time-invariant  to  an  observer  traveling  with  the 
crack  tip.  This  is  the  steady-state  case  when  depends  on  x  =  x^-a(t) 

only.  Now  we  have  that  a  =  constant  =  c  ,  x  =0,  and  (*)  =  -c„  d/dx. 

F  p  F 

The  solution  to  Eq.(37)  becomes 


dv,  .  0  B  C.  xz 

-  £  {-  £n(— )  -  -L  - 

dx  E  x  x  x  3 
P  XP 


(63) 
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Numerical  Results 

For  a  given  external  load,  we  presumably  know  in  terms  of  the 

crack  length  a(t).  A  relation  between  x  (t)  and  a(t)  can  subsequently 

P 

be  obtained  by  the  use  of  (57).  Hence,  in  principle,  a(t)  is  the  only 

unknown  quantity  in  (55).  An  equation  for  a(t)  and  a(t)  can,  for  example, 

be  obtained  from  (55)  by  the  use  of  the  critical  strain  criterion  for 

crack  propagation.  This  criterion  stipulates  that  crack  growth  will 

proceed  when  a  critical  strain  level  e  is  maintained  for  e  in  the 

cr  y 

plane  of  the  crack  at  a  characteristic  distance  A  ahead  of  the  crack  tip. 
It  appears,  however,  that  it  will  be  very  difficult  to  solve  a(t)  and 
a(t)  from  the  integral  equation  that  can  be  extracted  from  (55). 

In  the  examples  that  are  considered  here,  we  consider  an  inverse 
problem,  that  is,  we  prescribe  the  increasing  crack  length  a(t)  and 
the  variation  of  the  distant  tensile  stresses  o(t),  and  we  use  (55) 
and  (60)  to  compute  the  strain  at  a  small  distance  A  ahead  of  the 
crack  tip. 

Numerical  results  have  been  obtained  for  a  material  with  the 
following  mechanical  properties,  which  are  comparable  to  those  of 
CrMnSiNi  Steel: 

Young' 8  modulus:  E  =  2.06  x  1011  N/m2 

Poisson's  ratio:  v  =  0.3 

Yield  stress  in  shear:  k  ■  8.13  x  10®  N/m2 
Plane  stress  fracture  toughness:  K  *  16.7  x  107  N/m2 
The  geometry  considered  was  an  edge  crack  of  initial  length 
a 

o 


50mm 


in  a  half-plane.  The  half-plane  was  subjected  to  distant  tensile 
stresses  of  magnitude  a(t).  The  relevant  stress  intensity  factor  was 
taken  as 

Kj  =  1.1215  [xaCOl^aU)  0 

The  first  example  attempts  to  consider  a  case  where  the  stress 
intensity  factor  varies  in  such  a  manner  that  a  steady-state  situation 
can  be  established.  We  choose 


a(t)  =  ao 


t  <  t 
—  s 


0 


*  ao  +  a^t-tg)  +  a1(exp[-(t-tg) ]  -  1]}  t  >_  tg 

and 

o(t)  *  (t/t  )o  t  £  t 

O  U  O 


Here  oq  was  taken  as  the  time  that  K^,  reaches  the  value  of  the 
fracture  toughness 


0 

« 

0 

0 


o  -  K_/(1.1215/TTa  )  ,  0 

O  L  0 

while  t  is  the  time  that  the  crack  tip  arrives  at  a  distance  A  from 
c 

the  position  of  the  elastic-plastic  boundary  at  time  t  *  tg.  Thus 
t  can  be  computed  from 


a(t  )  +  A  *  a  +  x  (t  ) 


V  *  v*  h. "  O •  *  .*  .*  \  •-*  V  -V 


The  length  a^  was  taken  as  a^  =  63.9mm,  and  the  time  tg  was  taken 

as  tg  ■  30s.  It  follows  from  (67c)  and  (65)  that  remains  constant 

for  t  >  tc»  The  crack-tip  speed  a(t)  and  the  distant  tensile  stress 

o(t)  have  been  plotted  in  Fig.  2a  and  Fig.  2b,  respectively. 

The  strain  has  been  plotted  in  Fig.  3.  The  upper  curve 

represents  the  strain  at  a  fixed  position  A  =  1mm  ahead  of  the  crack 

tip.  For  t  <  t  =  30s,  this  is  a  fixed  position.  For  t  >  t  the 
s  s 

position  moves  with  the  crack  tip .  The  curves  numbered  2-10  represent 

the  strains  for  specific  fixed  material  points.  Points  2,  3  and  4 

were  located  inside  the  plastic  zone  at  time  t  »  t  ,  and  the  corres- 

s 

ponding  strains  were  computed  by  Eq.(60  with  upper  limit  t  <_  tg.  Points  5-10  were 

outside  the  plastic  zone  at  time  t  -  t  ,  and  for  these  points  the 

s 

strains  were  computed  by  Eq.(55).  Curves  2-10  all  end  at  a  time  t£ 
at  which  the  crack  tip  is  a  distance  A  from  the  material  point.  Time 
t£  is  computed  from  Eq.(61).  It  is  noted  that  the  strain  quickly 
approaches  an  apparent  steady-state  value,  which  is  just  equal  to  the 
steady-state  value  that  can  be  computed  from  Eq.(64).  Thus,  for  the 
distant  stress  given  by  Eq.(67),  stable  crack  propagation  at  a 
constant  strain  ey  is  established  quickly. 

For  the  second  example  we  choose 


o(t)  -  (7— )o  for  t  <  t 
to  —  s 

s 

■  o  for  t  >  t  , 
o  —  s 

where  is  defined  by  Eq.(68).  The  crack  starts  to  propagate  with 


(70a) 


(70b) 


(■  *>'•  .*•  .*•  .*•  ,%  .*•  .  .  V. S'  *.*  •. 
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constant  velocity  c^,  at  time  t  =  tg: 

a(t)  =  a  for  t  <  t  (71a) 

o  —  s 

=  a  +  c_(t-t  )  for  t  >  t  (71b) 

o  F  s  —  s 

For  three  crack-tip  speeds,  the  strains  at  a  location  which  remains  at 

a  fixed  distance  of  A  =  3mm  ahead  of  the  crack  tip  are  shown  in  Fig.  4. 

Thus,  the  location  where  the  strain  is  computed  is  stationary  for 

t  <  t  ,  and  it  moves  with  the  crack-tip  speed  c_  for  t  >  t  .  The 
—  s  r  s 

strains  increase  with  time,  and  the  rate  of  change  depends  significantly 
on  the  crack  tip  speed. 
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Fig.  4:  Strain  e  at  fixed  distance  A  ahead  of  crack  tip,  for  ramp  stress  defined 
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